Volume 99» Number 4, July-Augost J 994 " 

Journal of Research of the National Institute of Standards and Technology 

[J. Ras. Natl, Inst. Stand. Technol 9% 505 (1994)] 

Identification of Failure Origin Through 
Testing and the Weibull Risk-of-Rupture 



Volume 99 



Number 4 



July-^August 1994 



John A. Tesk, Martin Y. M. 
Chiang, and Spurgeon M, Keeny, 

m 

National Institute of Standards 
and Technology, 
Gaithersburg, MD 20899 

Jun Tang 

Shanghai Pharmaceutical 
Industry Design Institute, 
Shanghai, China 

and 

Yuuji Sato 

Hiroshima University School of 

Dentistry, 

Hiroshima, Japan 



The strcfis distribution in bond layers 
of two different tbickncsses (50 p.m and 
200 \Lm) was calculated by finite cle- 
ment anoly^is for pairs of rectangular 
erns.^ section metnl bars bonded to eaeh 
other and subjected tn four ptnni bend- 
ing. These stresses were used to aid \r\ 
identification of the fsiilury origin by 
use of the Weibiill tisk-of^rupiure (RR) 
function. By use of the stress distribu- 
tions, ihe characteristic strength from 
50 M-iTi l>ond test specimens eould be 
correlated with (bal for 200 jiin bond 
test j^pccimens when the failure was as- 
sumed to have an iiitorfacial origin. 
The finite clement meshes were refined 
twice and the ratios of characteristic 
strengths were recalculated and re- 
mained virtually unchanged l:y each of 
the mesh refincmcnls. tlcncc, tiie iden- 
tification of the intcrnite as the failure 



origin remained consistent. Further, the 
use of sitresscsi extrapolated to zero 
mesh iJize also produced the same ra- 
tios. Therefore, the RR calculations do 
not appear to be sensitive to the mesh 
sizes used for the stress calculations 
when the meshes are comparable or 
when changed in a compaf^ble manner. 
The results show this method can be 
consistent and a useful adjunct for 
idcntiflcatitm of failure origins. 
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1. Introduction 



In previous work [1,2] bending tests were con- 
ducted on adhesively bonded specimens of a dental 
alloy. The purpose was to determine; 

a) how much the bond thickness influenced the 
test results; 

b) whether the failure origins appeared to be the 
same for the two different bond thicknesses em- 
ployed; 

c) failure origin through analysis using the risk-of- 
rupture (RR) function. 



In this paper, the finite element method was used 
to arrive at the stress distributions used for the RR 
analyses that employ the Weibull risk-of- rupture 
function' [3]. An ancillary purpose, therefore^ was 
to ascertain how sensitive the analyses were to the 
fineness of the finite-element mesh and, hence, 
whether the method can be applied with confi- 
dence to the analyses conducted for a, b, and c 
above. 



'Today, a risk-of-rupture funeiiim, defined by Weibull, wtjuld 
be recognized as a hazard function. 
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2. Materials and Methods 

The failure of brittle materials is typically 
catastrophic and in many instances the failure 
stresses obtained from a set of test specimens fol- 
low a Weibull distribution. For a homogeneous 
isotropic material subjected to a uniform tensile 
stress, (T, the probability of failure, P((t), is given 
by: 



/>(o-) = l-exp-[5(trMur], 



(1) 



where aim is a characteristic strength for a speci- 
men of unit dimension; m is a shape factor 
{Weibull modulus); and 5 is a size factor (the ratio 
of the failure originating dimension to a unit di- 
mension of the same kind) and represents the vol- 
ume, K, area, A^ or other dimension in which 
reside the flaws from which the failure originates 
[3]. Eq. (1) is often written as: 



P(£r)=l^exp-[(cr/oa)'"L 



(la) 



where the size of the specimen, 8, is subsumed into 
oil. This form of the equation is commonly used 
when analyzing test data and the effects of speci- 
men size are ignored. It is clear from Eq, (1) that 
for specimens of two sizes, 8[ and §2, Viith the same 
failure origins and presenting the same distribution 
of failure stresses {mt=m2)^ there will be two dif- 
ferent values of oo for Eq, (la), with the larger size, 
call it 5i, leading to a value, ob^ that is less than 

For such specimen sets, the relation between the 
characteristic strengths calculated by Eq. (la) is 
[4]: 



*7b,2 = crij,i[5i/§2] 



Mm 



(2) 



For a nonuniform tensile stress field, a more gen- 
eral form of Eq. (1) is necessary: 



Pia)==l-tKp-UiitT/<r,)'"dS], 



(3) 



where the region of integration over 8 is the region 
critical to failure (rupture) and it can be in one* 
two, or three dimensions. Then the relationship be- 
tween the values of ^t^ [Eq. (la)] as determined 
from experiments on sets of specimens having ei- 
ther one or the other of the bond thicknesses, is 
given by the ratios of the exponents of Eq. (3), i.e., 






(4) 



d52 



Canceling terms in i7u.u on the right side and tron 
the left side leads to 






f (oi)""d52-' 



(4a) 



When o;(5) is not known as an explicit function, 
the relationship between au.i and a^j. may, in princi- 
ple, be approximated by [3]): 

ab.:-a-.u[(2^iA5i,)/(2:^iA52.0]^^ (5) 

and the validity of the approximation must be 
checked by compulation. Here the summations are 
over all the elements considered to be involved 
with the failure (interface, volume, etc.) and the 
stresses can be evaluated by the finite-element 
method of analysis- 
Note; The stress field in the bond region is typi- 
cally three dimensional; the analyses of this paper 
utilize unidirectional tensile stresses because alter- 
ations 10 the principal stresses were found to be 
minor and may be ignored. We also note that for a 
variety of reasons (plasticity, change in composi- 
tion, properties, or flaw populations) this analysis 
method may not apply for very thin bond layers 
approaching micrometers or less. 

Each assumed failure origin for a specimen has 
its own specific 5 with its associated stresses. When 
the ratios of volumes, surface areas, interface ar- 
eas, edge lengths etc, (any dimensions containing 
the flaws from which failures may originate) are 
properly chosen to be different for experimental 
tests, only one set of ft's, a/s, §'s, and oj's should 
produce coincidence between the experimentally 
determined ratio of oo's and the ratios of either the 
integrals shown by Eq, (4a) or the summations as 
shown in Eq. (5). 

As with any analysis employing the finite ele- 
ment (FE) method for determination of the 
stresses, a critical question arises as to the FE- 
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mesh sensitivity of Eq. (4). If sensitive, then the 
method would not, in actuality, be useful for the 
correlation of results from differently sized speci- 
mens. 

To provide insight into the ability of this ap- 
proach to identify sources of failure, rectangular 
bond specimens as shown in Fig. 1 were prepared 
for testing in four-|X)int bending, with either 50 iim 
or 200 |xm bond-thicknesses. Rectangular bars of 
the bulk bonding material were also tested in 
three-point bending. The details of specimen 
preparation were given in a presentation by Keeny 
et ai. [1]. The number of specimens and the results 
for each test series are shown in Table 1. 

A three-dimensional, finite-elemenl elasticity 
modeP was used for evaluation of the stress distri- 
bution throughout the volume of the bond region. 



The original bond model (Fig. 2) consisted of 2,197 
elements in one quadrant of the specimen which 
had three planes of symmetry. Subsequent refine* 
ments of interface and surface elements led to ele- 
ments 1/2 and 1/3 the original size. The validity of 
the model was checked by comparison of the finite- 
element results for a homogeneous beam with the 
analytical solution. Examples of how the bending 
tensile stresses change as a function of the thick- 
ness of the bond layer are shown in Figs. 3, 4, and 
5, for which E^=E^IS^ where: E\, is Young's mod- 
ulus for the bonding material; and £, is Young's 
modulus for the alloy. 

If the failure stresses are referenced to the 
stresses along the surface, the ratio of the operative 
(effective) dimensions, 5, from which the failures 
originate are given by Eq. (6): 




Fig. 1* Specimen used for bond testing. The shaded area repre- 
sents the bonding material between rectangular beams of alloy 
that were bonded together. Two small projections were used to 
control the width of the bond at cither 50 ^m or 2CX) \Lm. 



^itfj \ C7H, 



:)■ 



_J 



(6) 



where o-r denotes the reference stress. 

For these calculations the bending stress at the 
surface was used as the reference stress for calcula- 
tions of bending strength and (TRa-^A- Then, 
from Eq. (2) and Eq. (6) 



fT,Jo^j = (^cfu/5.fr,2)'"" = ^ ' 



Urn 



(7) 



which is equivalent to Eq. (5). 



Tabtc 1. Uncorrected Weibull parameters and mean strengths 



Test 



Gap (tim) N* cm (MPa) {range 95%f m {ran^e 95%} 



4-Pt Bend 



3-Pt Bend 



50 

200 

Bars of 
bulk bond 
material 



25 

25 

54 



110 {106-115} 
107 {103-111} 

SS {81-89} 



11.3 {8.2-14.0} 
11.518.3-14.2} 

(5.8 {5.5-7.93} 



*N = i\ic number of specimens. 

*'{} = the associated confidence bounds on vo and m as deit^rmined from the data. 



^Developed at NIST. 
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Adhesive 



Adhenend 



Fntcrfacc 



Ftg, 2* A 1/8 section of a three-dimensional mixlel for finite element calcula- 
tions of stress. The specimens (F[g. 1) had three p[ancs of s>Tnmctry, 



1.00 



0.90 



lo.eo 



0,70 



0,60 




50 100 1?0^ 200 

THICKNESS (MICROMETERS) 



Fig, 3. Result of finite ctcment calculation of the near-surface 
bending stress at the bond midplane that bisects the bend speci- 
mens into symmetrical halves. 



3, Results 

By the use of the right-hand side of Eq. (7) and 
the finite-element-derived stress distributions^ 
characteristic strength ratios were calculated for 
four potential regions (Table 2) where the failure 
of the bond could originate^ i.e., volume, surface, 
interface, and inierface-line-junction failure 
origins. These were then compared with the results 
obtained from the left hand side of Eq. (7). For 
these calculations, a value of w = U.5 was used 






Eo 



i.9 ■ 



0.6 




200 fM 



0.0 



S J I T1 1 1 rrT"P r t'r~t 



O.S 0.4 0.6 0.8 1.0 

NORMAUZED THICKNESS 



Fig. 4, Surface bending stress across ihe surface of the bond 
material, from one interface junaion to the opposite one, as 
calculated for ihrcc thicknesses. The 5 p.m thickness is pre- 
sented to illustrate the trend toward beam stress calculations as 
the thickness approaches zero. The deviations from beam thc- 
oiy calculations arc appreciable for thick bond specimens* show- 
ing the need to use the more robust finite clement method for 
the faihirc imulysis employing the RR function. 

[in Tabic 1, m was obtained from Newton-Raphson 
iteration for fitting experimental data to Eq. (1)]. 

The ratio of the experimental characteristic 
strengths has a 95% confidence range of 0.955 to 
LI L When this ratio is compared with ratios calcu- 
lated from the finite element analyses, the inter- 
face (Table 2) is identified as the origin of failures. 
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20 40 60 BO 100 120 140 160 160 200 
THICKNESS (MICROMETERS) 

Rg* 5- Finite element calculated stress for the adherend-adher' 
cnt interface surface junction line, illustrating the dramatie ef- 
fect of bond thickness on interface stresses. 



Table 1, Ratios of characteristic strengths calculated from risks 
nf rupture; (strength, 200 tim)/(sircngth, 50 p-m). Four-Point- 
bcndlng 



Assumed 
failure 

origin 


Coarse 
mesh 


Refined 
mesh 1 


Refined 
mesh 2 


Volume 


1.12 


1.12 


1.12 


Interface 


.995 


),0I 


L06 


Surface 


U26 


1.36 


1.32 


Interface- 
surface 
junction line 


1.24 


1.2S 


L41 


Interface- 
surface 
junction line 
(extrapolated)' 


138 


1.37 


1J7 



"Obtained from extrapolations of finite clement derived 
stresses at ccntroids to ihe interface between alloy and bond 
layer. 

with all other failure origins excluded. The row 
with the next closest match of strength ratios is that 
for volume failures which^ with a ratio of 1.12^ lies 
just outside the 95% confidence range, so this ar- 
gument, by itself, is somewhat unconvincing. How- 
ever, volume failures are ruled out because the m 
value of 11.5 from the bond tests differs, at the 
90% confidence level, from the value of m =6.8 
which was obtained from the bulk specimen test 



data. The bulk specimens can fail only by volume 
or surface failures* The strength ratio calculations 
for bond specimens rule out surface failures. The 
m value differences then are used to rule out vol- 
ume failures. 

Hence, the most reasonable explanation is that 
the bond specimens fail by interfacial failures. This 
is consistent with features of the failed specimens 
which always presented regions showing interfacial 
debonding. 

There is some possibility that the strength of the 
bond itself would depend on the bond thickness 
due to a change in material response (i.e., forma- 
tion of plasticity). Such effects obviously cannot be 
dealt with by the linear elastic analysis presented 
and within the confines of this analysis, interfa- 
cially initiated failure is concluded. 



4. Summarj' 

An analysis by Weibull RR for bonded speci- 
mens of two different sizes tested in bending has 

shown: 

1) Correlations between characteristic strengths, 
cTo's, were possible through the use of finite-ele- 
ment derived stresses in the RR analysis. 

2) The correlations were not sensitive to the par- 
ticular mesh chosen. 

3) For the interface, surface, and interface-surface 
junction line» the stress ratio calculations em- 
ploying the element centroid stresses are not 
significantly changed by use of stresses from ex- 
trapolations to the interface. The largest differ- 
ence is for the interface-surface junction line 
and these are shown in Table 2. 

4) Because the absolute magnitude of each RR 
calculation changes, mesh of the same size and 
configuration must be used for each set of com- 
parisons. 

5) The origins of failures can be determined by 
suitable testing and analysis of different size 
bond specimens and bulk specimens of the 
bonding material. This involves the use of a 
combined approach, analysis of the ao's and nt 
values. 

6) The determinalion of failure origin by this ap- 
proach can be useful for focusing attention on 
the proper parameters if improvements in sys- 
tem strength or performance are sought. 
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